Comparative Genomics Identifies Putative Interspecies Mechanisms Underlying Crbn-Sall4-Linked Thalidomide Embryopathy.

The identification of thalidomide-Cereblon-induced SALL4 degradation has brought new understanding for thalidomide embryopathy (TE) differences across species. Some questions, however, regarding species variability, still remain. The aim of this study was to detect sequence divergences between species, affected or not by TE, and to evaluate the regulated gene co-expression in a murine model. Here, we performed a comparative analysis of proteins experimentally established as affected by thalidomide exposure, evaluating 14 species. The comparative analysis, regarding synteny, neighborhood, and protein conservation, was performed in 42 selected genes. Differential co-expression analysis was performed, using a publicly available assay, GSE61306, which evaluated mouse embryonic stem cells (mESC) exposed to thalidomide. The comparative analyses evidenced 20 genes in the upstream neighborhood of NOS3, which are different between the species who develop, or not, the classic TE phenotype. Considering protein sequence alignments, RECQL4, SALL4, CDH5, KDR, and NOS2 proteins had the biggest number of variants reported in unaffected species. In co-expression analysis, Crbn was a gene identified as a driver of the co-expression of other genes implicated in genetic, non-teratogenic, limb reduction defects (LRD), such as Tbx5, Esco2, Recql4, and Sall4; Crbn and Sall4 were shown to have a moderate co-expression correlation, which is affected after thalidomide exposure. Hence, even though the classic TE phenotype is not identified in mice, a deregulatory Crbn-induced mechanism is suggested in this animal. Functional studies are necessary, especially evaluating the genes responsible for LRD syndromes and their interaction with thalidomide-Cereblon.

Methylphenidate effects on mice odontogenesis and connections with human odontogenesis.

The purpose of this study is to evaluate the effects of Methylphenidate exposure on mice odontogenesis and connect them by bioinformatics with human odontogenesis. Thirty-two pregnant Swiss mice were divided into treated group and control group, which received, respectively, 5 mg/kg of Methylphenidate and saline solution from the 5th to the 17th day of pregnancy. The mouse embryos tooth germs were analyzed through optical microscopy, and the data collected were analyzed statistically by Fisher's exact test. The presence and similarity of Methylphenidate-associated genes (Pharmgkb database) in both organisms and their interaction with dental development genes (AmiGO2 database) were verified on STRING database. Rates of tooth germ malformations were higher in treated than in control group (Control: 18; Treated: 27; p = 0.035). Mouse embryo malformations were connected with 238 interactions between 69 dental development genes with 35 Methylphenidate genes. Fourteen interactions for four Methylphenidate genes with four dental development genes, with human experimental data, were connected with mouse phenotype data. By homology, the interactions and conservation of proteins/genes may indicate similar outcomes for both organisms. The exposure to Methylphenidate during pregnancy affected odontogenesis in mouse embryos and may affect human odontogenesis. The study of malformations in mice, with a bioinformatics approach, could contribute to understanding of the Methylphenidate effect on embryo development. These results may provide novel hypotheses for further testing and reinforce the FDA protocol: as Methylphenidate is included in category C, its use during pregnancy should be considered if the benefits outweigh the risks.

Genetic variants in miRNAs differentially expressed during brain development and their relevance to psychiatric disorders susceptibility.

OBJECTIVES: MicroRNAs (miRNAs) play an important regulatory role in the expression of genes involved in brain functions during development. Genetic variants in miRNA genes may impact their regulatory function and lead to psychiatric disorders. To evaluate the role of genetic variants in genes of miRNAs differentially expressed during neurodevelopment on autism spectrum disorder (ASD), attention deficit/hyperactivity disorder (ADHD), schizophrenia (SCZ), and major depressive disorder (MDD). METHODS: The miRNAs were identified in the literature. Summary statistics from the most recent genome-wide association studies to date were used to evaluate the association between the selected polymorphisms and each disorder in a look-up approach. In a global analysis, we compared the standardised risk effect of variants in neurodevelopment-related miRNAs with those in the remaining miRNAs from miRBase. RESULTS: The global analysis showed that variants in neurodevelopment-related miRNAs had higher risk effects compared to the other miRNAs for SCZ (p = 0.010) and ADHD (p = 0.001). MIR33B, MIR29B2, MIR29C, MIR137, and MIR135A1 were significantly associated with SCZ, while 55.9% of the miRNAs were at least nominally associated with one or more psychiatric disorders (p < 0.05). CONCLUSIONS: Genetic variants in neurodevelopment-related miRNAs play an important role in the genetic susceptibility of psychiatric disorders, mainly SCZ and ADHD.

Evolutionary analysis of the Musashi family: What can it tell us about Zika?

Despite worldwide research efforts since 2015, Zika virus infection and its consequences are not fully understood yet. Nowadays, it is known that microcephaly is only one of the possible outcomes of being infected by ZIKV during the early stages of life. Musashi 1 (MSI1) is an RNA-binding protein that is involved in neurodevelopmental processes. Also, ZIKV genome (a single-stranded positive-sense RNA) uses MSI1 for its replication. Here we perform an evolutionary analysis of MSI1 coding sequence and their orthologs in vertebrate species. We added original sequencing data from selected regions of interest (RNA-binding domains-RBDs of MSI1) of sixteen Platyrrhini (or New World monkeys), known to have high evolutionary rates. The Musashi family (MF) includes MSI2, TARDBP, DAZAP1, HNRNPD, HNRNPDL, and HNRNPAB, which do not interact with the virus but are critical RNA-binding proteins that act on many regulatory processes ubiquitously. We found that all sixteen primate species have the RBD1 of MSI1 conserved. While the general code sequences of MF genes are under purifying selection, the evolution of regulatory mechanisms, especially alternative splicing, seems to be a frequent phenomenon in these genes. Different isoforms differ in the N-terminal region and it affects protein size. Existing MSI1 isoforms probably diverge in their binding affinity, the kinetics of interaction, and other aspects when in the MSI1-ZIKV complex. It is a signal that some RBD-containing MSI1 isoforms can be incompatible to ZIKV binding and replication. Consequently, the chance of ZIKV successfully infecting host cells could also be associated with alternative splicing and expression of ZIKV-compatible MSI1 isoforms in both inter and intraspecific levels.

CRL4-Cereblon complex in Thalidomide Embryopathy: a translational investigation.

The Cereblon-CRL4 complex has been studied predominantly with regards to thalidomide treatment of multiple myeloma. Nevertheless, the role of Cereblon-CRL4 in Thalidomide Embryopathy (TE) is still not understood. Not all embryos exposed to thalidomide develop TE, hence here we evaluate the role of the CRL4-Cereblon complex in TE variability and susceptibility. We sequenced CRBN, DDB1, CUL4A, IKZF1, and IKZF3 in individuals with TE. To better interpret the variants, we suggested a score and a heatmap comprising their regulatory effect. Differential gene expression after thalidomide exposure and conservation of the CRL4-Cereblon protein complex were accessed from public repositories. Results suggest a summation effect of Cereblon variants on pre-axial longitudinal limb anomalies, and heatmap scores identify the CUL4A variant rs138961957 as potentially having an effect on TE susceptibility. CRL4-Cereblon gene expression after thalidomide exposure and CLR4-Cereblon protein conservation does not explain the difference in Thalidomide sensitivity between species. In conclusion, we suggest that CRL4-Cereblon variants act through several regulatory mechanisms, which may influence CRL4-Cereblon complex assembly and its ability to bind thalidomide. Human genetic variability must be addressed not only to further understand the susceptibility to TE, but as a crucial element in therapeutics, including in the development of pharmacogenomics strategies.

Clusters of genetic diseases in Brazil.

The aim of this paper is to present a database of isolated communities (CENISO) with high prevalence of genetic disorders or congenital anomalies in Brazil. We used two strategies to identify such communities: (1) a systematic literature review and (2) a "rumor strategy" based on anecdotal accounts. All rumors and reports were validated in a stepwise process. The bibliographical search identified 34 rumors and 245 rumors through the rumor strategy, and 144 were confirmed. A database like this one presented here represents an important tool for the planning of health priorities for rare diseases in low- and middle-income countries with large populations.

Searching for ancient balanced polymorphisms shared between Neanderthals and Modern Humans.

Hominin evolution is characterized by adaptive solutions often rooted in behavioral and cognitive changes. If balancing selection had an important and long-lasting impact on the evolution of these traits, it can be hypothesized that genes associated with them should carry an excess of shared polymorphisms (trans- SNPs) across recent Homo species. In this study, we investigate the role of balancing selection in human evolution using available exomes from modern (Homo sapiens) and archaic humans (H. neanderthalensis and Denisovan) for an excess of trans-SNP in two gene sets: one associated with the immune system (IMMS) and another one with behavioral system (BEHS). We identified a significant excess of trans-SNPs in IMMS (N=547), of which six of these located within genes previously associated with schizophrenia. No excess of trans-SNPs was found in BEHS, but five genes in this system harbor potential signals for balancing selection and are associated with psychiatric or neurodevelopmental disorders. Our approach evidenced recent Homo trans-SNPs that have been previously implicated in psychiatric diseases such as schizophrenia, suggesting that a genetic repertoire common to the immune and behavioral systems could have been maintained by balancing selection starting before the split between archaic and modern humans.

Searching for ancient balanced polymorphisms shared between Neanderthals and Modern Humans

Abstract Hominin evolution is characterized by adaptive solutions often rooted in behavioral and cognitive changes. If balancing selection had an important and long-lasting impact on the evolution of these traits, it can be hypothesized that genes associated with them should carry an excess of shared polymorphisms (trans- SNPs) across recent Homo species. In this study, we investigate the role of balancing selection in human evolution using available exomes from modern (Homo sapiens) and archaic humans (H. neanderthalensis and Denisovan) for an excess of trans-SNP in two gene sets: one associated with the immune system (IMMS) and another one with behavioral system (BEHS). We identified a significant excess of trans-SNPs in IMMS (N=547), of which six of these located within genes previously associated with schizophrenia. No excess of trans-SNPs was found in BEHS, but five genes in this system harbor potential signals for balancing selection and are associated with psychiatric or neurodevelopmental disorders. Our approach evidenced recent Homo trans-SNPs that have been previously implicated in psychiatric diseases such as schizophrenia, suggesting that a genetic repertoire common to the immune and behavioral systems could have been maintained by balancing selection starting before the split between archaic and modern humans.

Functional New World monkey oxytocin forms elicit an altered signaling profile and promotes parental care in rats.

The neurohormone oxytocin is a key player in the modulation of reproductive and social behavioral traits, such as parental care. Recently, a correlation between different forms of oxytocin and behavioral phenotypes has been described in the New World Monkeys (NWMs). Here, we demonstrate that, compared with the Leu8OXT found in most placental mammals, the Cebidae Pro8OXT and Saguinus Val3Pro8OXT taxon-specific variants act as equi-efficacious agonists for the Gq-dependent pathway but are weaker agonists for the ß-arrestin engagement and subsequent endocytosis toward the oxytocin receptor (OXTR). Upon interaction with the AVPR1a, Pro8OXT and the common Leu8OXT yielded similar signaling profiles, being equally efficacious on Gq and ß-arrestin, while Val3Pro8OXT showed reduced relative efficacy toward ß-arrestin. Intranasal treatment with either of the variants increased maternal behavior and also promoted unusual paternal care in rats, as measured by pup-retrieval tests. We therefore suggest that Val3Pro8OXT and Pro8OXT are functional variants, which might have been evolutionarily co-opted as an essential part of the adaptive genetic repertoire that allowed the emergence of taxon-specific complex social behaviors, such as intense parental care in the Cebidae and the genus Saguinus.

FOXP in Tetrapoda: Intrinsically Disordered Regions, Short Linear Motifs and their evolutionary significance.

The FOXP subfamily is probably the most extensively characterized subfamily of the forkhead superfamily, playing important roles in development and homeostasis in vertebrates. Intrinsically disorder protein regions (IDRs) are protein segments that exhibit multiple physical interactions and play critical roles in various biological processes, including regulation and signaling. IDRs in proteins may play an important role in the evolvability of genetic systems. In this study, we analyzed 77 orthologous FOXP genes/proteins from Tetrapoda, regarding protein disorder content and evolutionary rate. We also predicted the number and type of short linear motifs (SLIMs) in the IDRs. Similar levels of protein disorder (approximately 70%) were found for FOXP1, FOXP2, and FOXP4. However, for FOXP3, which is shorter in length and has a more specific function, the disordered content was lower (30%). Mammals showed higher protein disorders for FOXP1 and FOXP4 than non-mammals. Specific analyses related to linear motifs in the four genes showed also a clear differentiation between FOXPs in mammals and non-mammals. We predicted for the first time the role of IDRs and SLIMs in the FOXP gene family associated with possible adaptive novelties within Tetrapoda. For instance, we found gain and loss of important phosphorylation sites in the Homo sapiens FOXP2 IDR regions, with possible implication for the evolution of human speech.

FOXP in Tetrapoda: Intrinsically Disordered Regions, Short Linear Motifs and their evolutionary significance

Abstract The FOXP subfamily is probably the most extensively characterized subfamily of the forkhead superfamily, playing important roles in development and homeostasis in vertebrates. Intrinsically disorder protein regions (IDRs) are protein segments that exhibit multiple physical interactions and play critical roles in various biological processes, including regulation and signaling. IDRs in proteins may play an important role in the evolvability of genetic systems. In this study, we analyzed 77 orthologous FOXP genes/proteins from Tetrapoda, regarding protein disorder content and evolutionary rate. We also predicted the number and type of short linear motifs (SLIMs) in the IDRs. Similar levels of protein disorder (approximately 70%) were found for FOXP1, FOXP2, and FOXP4. However, for FOXP3, which is shorter in length and has a more specific function, the disordered content was lower (30%). Mammals showed higher protein disorders for FOXP1 and FOXP4 than non-mammals. Specific analyses related to linear motifs in the four genes showed also a clear differentiation between FOXPs in mammals and non-mammals. We predicted for the first time the role of IDRs and SLIMs in the FOXP gene family associated with possible adaptive novelties within Tetrapoda. For instance, we found gain and loss of important phosphorylation sites in the Homo sapiens FOXP2 IDR regions, with possible implication for the evolution of human speech.

The role of protein intrinsic disorder in major psychiatric disorders.

Although new candidate genes for Autism Spectrum Disorder (ASD), Schizophrenia (SCZ), Attention-Deficit/Hyperactivity Disorder (ADHD), and Bipolar Disorder (BD) emerged from genome-wide association studies (GWAS), their underlying molecular mechanisms remain poorly understood. Evidences of the involvement of intrinsically disordered proteins in diseases have grown in the last decade. These proteins lack tridimensional structure under physiological conditions and are involved in important cellular functions such as signaling, recognition and regulation. The aim of the present study was to identify the role and abundance of intrinsically disordered proteins in a set of psychiatric diseases and to test whether diseases are different regarding protein intrinsic disorder. Our hypothesis is that differences across psychiatric illnesses phenotypes and symptoms may arise from differences in intrinsic protein disorder content and properties of each group. A bioinformatics prediction of intrinsic disorder was performed in proteins retrieved based on top findings from GWAS, Copy Number Variation and candidate gene investigations for each disease. This approach revealed that about 80% of studied proteins presented long stretches of disorder. This amount was significantly higher than that observed in general eukaryotic proteins, and those involved in cardiovascular diseases. These results suggest that proteins with intrinsic disorder are a common feature of neurodevelopment and synaptic transmission processes which are potentially involved in the etiology of psychiatric diseases. Moreover, we identified differences between ADHD and ASD when the binary prediction of structure and putative binding sites were compared. These differences may be related to variation in symptom complexity between both diseases. © 2016 Wiley Periodicals, Inc.

The contribution of protein intrinsic disorder to understand the role of genetic variants uncovered by autism spectrum disorders exome studies.

Several autism spectrum disorders (ASD) exome studies suggest that coding single nucleotide variants (SNVs) play an important role on ASD etiology. Usually, the pathogenic effect of missense mutations is estimated through predictors that lose accuracy for those SNVs placed in intrinsically disordered regions of protein. Here, we used bioinformatics tools to investigate the effect of mutations described in ASD published exome studies (549 mutations) in protein disorder, considering post-translational modification, PEST and Molecular Recognition Features (MoRFs) motifs. Schizophrenia and type 2 diabetes (T2D) datasets were created for comparison purposes. The frequency of mutations predicted as disordered was comparable among the three datasets (38.1% in ASD, 35.7% in schizophrenia, 46.4% in T2D). However, the frequency of SNVs predicted to lead a gain or loss of functional sites or change intrinsic disorder tendencies was higher in ASD and schizophrenia than T2D (46.9%, 36.4%, and 23.1%, respectively). The results obtained by SIFT and PolyPhen-2 indicated that 38.9% and 34.4% of the mutations predicted, respectively, as tolerated and benign showed functional alterations in disorder properties. Given the frequency of mutations placed in IDRs and their functional impact, this study suggests that alterations in intrinsic disorder properties might play a role in ASD and schizophrenia etiologies. They should be taken into consideration when researching the pathogenicity of mutations in neurodevelopmental and psychiatric diseases. Finally, mutations with functional alterations in disorder properties must be potential targets for in vitro and in vivo functional studies.

Evolutionary pattern in the OXT-OXTR system in primates: coevolution and positive selection footprints.

Oxytocin is a nonapeptide involved in a wide range of physiologic and behavioral functions. Until recently, it was believed that an unmodified oxytocin sequence was present in all placental mammals. This study analyzed oxytocin (OXT) in 29 primate species and the oxytocin receptor (OXTR) in 21 of these species. We report here three novel OXT forms in the New World monkeys, as well as a more extensive distribution of a previously described variant (Leu8Pro). In structural terms, these OXTs share the same three low-energy conformations in solution during molecular dynamic simulations, with subtle differences in their side chains. A consistent signal of positive selection was detected in the Cebidae family, and OXT position 8 showed a statistically significant (P = 0.013) correlation with litter size. Several OXTR changes were identified, some of them promoting gain or loss of putative phosphorylation sites, with possible consequences for receptor internalization and desensitization. OXTR amino acid sites are under positive selection, and intramolecular and intermolecular coevolutionary processes with OXT were also detected. We suggest that some New World monkey OXT-OXTR forms can be correlated to male parental care through the increase of cross-reactivity with its correlated vasopressin system.

Evolutionary history of chordate PAX genes: dynamics of change in a complex gene family.

Paired box (PAX) genes are transcription factors that play important roles in embryonic development. Although the PAX gene family occurs in animals only, it is widely distributed. Among the vertebrates, its 9 genes appear to be the product of complete duplication of an original set of 4 genes, followed by an additional partial duplication. Although some studies of PAX genes have been conducted, no comprehensive survey of these genes across the entire taxonomic unit has yet been attempted. In this study, we conducted a detailed comparison of PAX sequences from 188 chordates, which revealed restricted variation. The absence of PAX4 and PAX8 among some species of reptiles and birds was notable; however, all 9 genes were present in all 74 mammalian genomes investigated. A search for signatures of selection indicated that all genes are subject to purifying selection, with a possible constraint relaxation in PAX4, PAX7, and PAX8. This result indicates asymmetric evolution of PAX family genes, which can be associated with the emergence of adaptive novelties in the chordate evolutionary trajectory.

The role of a mineralocorticoid receptor gene functional polymorphism in the symptom dimensions of persistent ADHD.

Attention-deficit/hyperactivity disorder (ADHD) affects approximately 5 % of school-aged children and 2.5 % of adults. Genetic studies in ADHD have pointed to genes in different neurobiological systems, with relatively small individual effects. The mineralocorticoid receptor is the main receptor involved in the initial triggering of stress response. Therefore, its encoding gene (NR3C2) is a candidate for psychiatric disorder studies, including ADHD, and behavioral phenotypes. There is evidence that the Val allele of the MRI180V polymorphism (rs5522) increases the risk of depression, attention and cognitive deficits. We investigated the possible role of the mineralocorticoid receptor gene in the symptom dimensions and susceptibility to persistent ADHD. We compared genotype and allele frequencies in 478 adult patients with ADHD and 597 controls and symptom dimensions in 449 patients and 132 controls. Diagnoses were based on the DSM-IV criteria. ADHD symptom dimensions were investigated with SNAP-IV for ADHD severity and Barkley scales for severity and impairment. Carriers of the Val allele presented higher inattention, hyperactivity/impulsivity and impairment scores, while genotype and allele frequencies did not differ between patients and controls. These results are consistent with a possible link between genetic variations in the HPA axis and inattention and hyperactivity measures.

Evolutionary responses to a constructed niche: ancient Mesoamericans as a model of gene-culture coevolution.

Culture and genetics rely on two distinct but not isolated transmission systems. Cultural processes may change the human selective environment and thereby affect which individuals survive and reproduce. Here, we evaluated whether the modes of subsistence in Native American populations and the frequencies of the ABCA1*Arg230Cys polymorphism were correlated. Further, we examined whether the evolutionary consequences of the agriculturally constructed niche in Mesoamerica could be considered as a gene-culture coevolution model. For this purpose, we genotyped 229 individuals affiliated with 19 Native American populations and added data for 41 other Native American groups (n = 1905) to the analysis. In combination with the SNP cluster of a neutral region, this dataset was then used to unravel the scenario involved in 230Cys evolutionary history. The estimated age of 230Cys is compatible with its origin occurring in the American continent. The correlation of its frequencies with the archeological data on Zea pollen in Mesoamerica/Central America, the neutral coalescent simulations, and the F(ST)-based natural selection analysis suggest that maize domestication was the driving force in the increase in the frequencies of 230Cys in this region. These results may represent the first example of a gene-culture coevolution involving an autochthonous American allele.

Glutamate and synaptic plasticity systems and smoking behavior: results from a genetic association study.

Smoking behavior is a multifactorial phenotype with significant heritability. Identifying the specific loci that influence smoking behavior could provide important etiological insights and facilitate the development of treatments to further reduce smoking related mortality. Although several studies pointed to different candidate genes for smoking, there is still a need for replication especially in samples from different countries. In the present study, we investigated whether 21 positive signals for smoking behavior from these studies are replicated in a sample of 531 blood donors from the Brazilian population. The polymorphisms were chosen based on their representativeness of different candidate biologic systems, strength of previous evidence, location and allele frequencies. By genotyping with the Sequenom MassARRAY iPLEX platform and subsequent statistical analysis using Plink software, we show that two of the SNPs studied, in the SLC1A2 (rs1083658) and ACTN1 (rs2268983) genes, were associated with smoking behavior in our study population. These genes are involved in crucial aspects of nicotine dependence, glutamate system and synaptic plasticity, and as such, are biologically plausible candidates that merit further molecular analyses so as to clarify their potential role in smoking behavior.

Linking dopamine neurotransmission and neurogenesis: The evolutionary history of the NTAD (NCAM1-TTC12-ANKK1-DRD2) gene cluster.

Genetic studies have long suggested the important role of the DRD2 gene in psychiatric disorders and behavior. Further research has shown a conjoined effect of genes in the Chr11q22-23 region, which includes the NCAM1, TTC12, ANKK1 and DRD2 genes, or NTAD cluster. Despite a growing need to unravel the role of this cluster, few studies have taken into account interspecies and evolutionary approaches. This study shows that behaviorally relevant SNPs from the NTAD cluster, such as rs1800497 (Taq1A) and rs6277, are ancient polymorphisms that date back to the common ancestor between modern humans and Neanderthals/Denisovans. Conserved synteny and neighborhood indicate the NTAD cluster seems to have been established at least 400 million years ago, when the first Sarcopterygians emerged. The NTAD genes are apparently co-regulated and this could be attributed to adaptive functional properties, including those that emerged when the central nervous system became more complex. Finally, our findings indicate that NTAD genes, which are related to neurogenesis and dopaminergic neurotransmission, should be approached as a unit in behavioral and psychiatric genetic studies.

The TP53 fertility network.

The TP53 gene, first described in 1979, was identified as a tumor suppressor gene in 1989, when it became clear that its product, the p53 nuclear phosphoprotein, was frequently inactivated in many different forms of cancers. Nicknamed "guardian of the genome", TP53 occupies a central node in stress response networks. The p53 protein has a key role as transcription factor in limiting oncogenesis through several growth suppressive functions, such as initiating apoptosis, senescence, or cell cycle arrest. The p53 protein is directly inactivated in about 50% of all tumors as a result of somatic gene mutations or deletions, and over 80% of tumors demonstrate dysfunctional p53 signaling. Beyond the undeniable importance of p53 as a tumor suppressor, an increasing number of new functions for p53 have been reported, including its ability to regulate energy metabolism, to control autophagy, and to participate in various aspects of differentiation and development. Recently, studies on genetic variations in TP53 among different populations have led to the notion that the p53 protein might play an important role in regulating fertility. This review summarizes current knowledge on the basic functions of different genes of the TP53 family and TP53 pathway with respect to fertility. We also provide original analyses based on genomic and genotype databases, providing further insights into the possible roles of the TP53 pathway in human reproduction.